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A fully first principles theory capable of treating strongly correlated solids remains the outstanding 
challenge of modern day materials science. This is exemplified by the transition metal oxides, 
prototypical Mott insulators, that remain insulating even in the absence of long range magnetic 
order. Capturing the non-magnetic insulating state of these materials presents a difficult challenge 
for any modern electronic structure theory. In this paper we demonstrate that reduced density 
matrix functional theory, in conjunction with the power functional, can successfully treat the non¬ 
magnetic insulating state of the transition metal oxides NiO and MnO. We show that the electronic 
spectrum retains a gap even in the absence of spin order. We further discuss the detailed way 
in which RDMFT performs for Mott insnlators and band insnlators, finding that for the latter 
occnpation number minimization alone is reqnired, but for the former full minimization over both 
occnpation numbers and natural orbitals is essential. 

PACS numbers: 71.10.-w, 71.27.+a, 71.45.Gm, 71.20.Nr 


INTRODUCTION 

One of the most useful derivatives of a ground- 
state density functional theory (DFT) calculation is 
the Kohn-Sham (KS) eigenvalues [T], which lead to a 
non-interacting spectrum. Even though the KS equa¬ 
tions represent an auxiliary non-interacting system whose 
states and eigenvalues may be quite different from the 
true quasi-particle system, empirical evidence shows 
that in many cases this single particle KS spectrum 
is in agreement with the x-ray photo-emission Spec¬ 
troscopy (XPS) and Bremsstrahlung isochromat spec¬ 
troscopy (BIS) experiments [5H5]. However, for strongly 
correlated materials this KS spectrum is found to be in 
fundamental disagreement with experimental reality. In 
the absence of spin-ordering all modern exchange cor¬ 
relation (xc) functionals within DFT fail to predict an 
insulating ground-state for transition metal mono-oxides 
(TMOs), the prototypical Mott insulators. On the other 
hand, it is well known experimentally that these materi¬ 
als are insulating in nature even at elevated temperatures 
(much above the Neel temperature) [6l [7] . This indicates 
that magnetic order is not the driving mechanism for the 
gap, but merely a co-occurring phenomenon. In fact not 
only DFT, but most modern many-body techniques such 
as the GW method also fail to capture the insulating be¬ 
havior in TMOs without explicit long range spin ordering 

[SHIQ]. 

In this regard, the two many-body techniques that are 
able to capture the correct physics of strong correlations 
are dynamical mean field theory (DMFT)[TTHl3] and re¬ 
duced density matrix functional theory (RDMFT) |14j: 
these two methods predicts TMOs as insulators, even in 


the absence of long range spin-order. This clearly points 
towards the ability of these techniques to capture physics 
well beyond the reach of most modern day ground-state 
methods. 


Despite this success the effectiveness of RDMFT as 
a ground state theory has been seriously hampered due 
to the absence of a technique for the determination of 
the spectral information. Recently, this final hurdle has 
also been removed and the spectral information thus ob¬ 
tained for TMOs was shown to be in good agreement with 
experiments [in]. However, these spectra were calculated 
in the presence of anti-ferromagnetic order. The question 
then arises as to how effective RDMFT is in describing 
the insulating state of Mott insulators in the absence of 
long range spin order. In order to answer this question, 
in the present work, we study the spectral properties of 
non-magnetic NiO and MnO. Here former is insulating 
due to interplay of Mott localization and charge transfer 
effects while the latter is insulating purely due to strong 
Mott localization. A detailed analysis of RDMFT and 
KS orbitals is performed which shows that, unlike in the 
case of band insulators, for Mott insulators the nature of 
two set of orbitals are very different and this difference 
is indeed crucial for the success of RDMFT in describing 
Mott physics. 
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THEORY 

Within RDMFT the one-body reduced density matrix 
( 1 -RDM) is the basic variable [TBIITT] 


All calculations are performed using the full-potential 
linearized augmented plane wave code Elk|40j. with prac¬ 
tical details of the calculations following the schemes de¬ 
scribed in Refs. M and HSl 


7 (r,r') = N 



.. drN^*ir', r2 ... rAr)$(r, r2 ... rjv), 

( 1 ) 


where $ denotes the many-body wave function. Di- 
agonalization of this matrix produces a set of natural 
orbitalsfTB]. (pj^, and occupation numbers, rij^, leading 
to the spectral representation 


7(1-, r') = (2) 

i.k 


where the necessary and sufficient conditions for ensem¬ 
ble iV-representability of 7 [15] require 0 < Uj^ < 1 for 
all j, k, and k = N. Here j represents the band 
index and k the crystal momentum. 

In terms of 7 , the total ground state energy m of the 
interacting system is (atomic units are used throughout) 
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where p{r) = 7 (r,r), Vext is a given external potential, 
and Axe we call the xc energy functional. In principle, 
Gilbert’s [T7| generalization of the Hohenberg-Kohn the¬ 
orem to the 1-RDM guarantees the existence of a func¬ 
tional A [ 7 ] whose minimum yields the exact 7 and the ex¬ 
act ground-state energy of systems characterized by the 
external potential I4xt(r). In practice, however, the cor¬ 
relation energy is an unknown functional of the 1-RDM 
and must be approximated. Although there are several 
known approximations for the xc energy functional [151- 
134] , the most promising for extended systems is the power 
functional [CT [TS] where the xc energy reads 

AxcH = i?xc[{<(>.k}, M] =-\J J 

(4) 


here 7 “ indicates the power used in the operator sense 

i.e. 


i 

for a = 1/2 this is the Muller functional [35], which is 
known to severely overestimate electron correlation [33]- 
139] while for a = 1 this functional is equivalent to the 
Hartree-Fock method, which includes no correlations. If 
a is chosen to be 1/2 < a < 1 , the power functional 
interpolates between the uncorrelated Hartree-Fock limit 
and the over-correlating Muller functional. 


RESULTS 



FIG. 1: (color online) Density of states as a function of energy 
(in eV) for NiO (left panel) and MnO (right panel). Results 
are obtained with (black) and without (red) long range (anti- 
ferromagnetic) spin order. For comparison experimental data 
taken from Refs. [3] and [5] is also shown (grey shaded area). 
Chemical potential is shown at dotted vertical line 

Presented in Fig. [T]are the spectra for the Mott insula¬ 
tors under consideration. It is immediately apparent that 
RDMFT captures the essence of Mott-Hubbard physics: 
both NiO and MnO present substantial gaps at the Fermi 
energy and are thus insulating in the absence of spin or¬ 
der. This fact was already noticed in a previous work [14] 
in which the presence of gap without any spin-order was 
deduced via very different technique, namely the discon¬ 
tinuity in the chemical potential as a function of the par¬ 
ticle number. A comparison of the non-magnetic spectra 
with the experimental data shows that the shape of the 
conduction band is well reproduced for both materials, 
but that the shape of the valence band is not in very 
good agreement with experiments. This agreement im¬ 
proves on inclusion of the spin order, indicating that even 
though insulating nature of TMO’s is not driven by spin 
order, spin polarization significantly effects the spectra 
of these materials. This is hardly surprising given that 
NiO and MnO have very large local moments of 1.9/iB 
and A.lpB respectively. 

Correct treatment of correlations is crucial for TMOs, 
the prototypical strongly correlated materials. As men¬ 
tioned above the power functional interpolates between 
two limits - the highly over correlated Muller [a = 0.5) 
and totally uncorrelated Hartree-Fock (a = 1). We now 
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FIG. 2: Density of states as a function of energy (in eV) 
for NiO (upper panels), MnO (middle panels), and Si (lower 
panels). The results are obtained using different values of a 
in Eq. 


look at the effect of correlations, by varying a, on the 
spectra of Mott insulators (NiO and MnO) and band 
insulator (Si), see Fig. The behaviour of the spec¬ 
tra as a function of a is rather trivial for band insula¬ 
tor, Si; the valence bands rigidly shift lower in energy 
leading to increase in the band gap. The behaviour for 
Mott insulators is different in that the shape of the bands 
change as a function of a. Both for NiO and MnO over 
correlated Muller functional incorrectly gives a metallic 
ground-state. For NiO, which has even number of elec¬ 
trons in a unit-cell, the Hartree-Fock method leads to a 
very large band gap insulator. In contrast to this, for 
MnO, with odd number of electrons in the unit-cell, a 
single particle theory such as Hartree-Fock can only give 
rise to a metallic ground state. This leads to highly non 
trivial behaviour for MnO as a function of a, which must 
lie within a small range (between 0.65 and 0.7) in which 
the correct insulating ground-state is obtained. Reassur¬ 
ingly, this is also the range of a in which correct ground 
state behaviour is seen for NiO. 



Energy (eV) 

FIG. 3: Density of states as a function of energy (in eV) for 
NiO (top panel), MnO (middle panel) and Si (lower panel). 
Results are obtained with (black) and without (red) opti¬ 
mization of the natural orbitals with in RDMFT. KS results 
(green) are obtained using local density approximation [41]. 

Within RDMFT there are no Kohn-Sham-like equa¬ 
tions to solve, and a direct minimization over natural 
orbitals and occupation numbers is required while main¬ 
taining the ensemble 7V-representability conditions. The 
minimization over occupation numbers is computation¬ 
ally very efficient (for details see Ref. |T4|), but the same 
cannot be said about the minimization over the natu¬ 
ral orbitals. In practical terms, the natural orbitals (see 
Eq. Q) are expanded in a set of previously converged 
KS states, and optimization of the natural orbitals is 
performed by varying the expansion coefficients. This 
procedure allows us to examine how different KS states 
are from fully optimized natural orbitals. In the present 
work these KS states were obtained using local density 
approximation (LDA)[1T|. 
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In Fig. three set of results are shown; (i) KS density 
of states, (ii) RDMFT density of states obtained without 
optimizing the natural orbitals i.e. by using KS orbitals 
as natural orbitals but fully optimizing the occupation 
numbers and (iii) the fully optimized RDMFT results 
i.e. full optimization over the natural orbitals and occu¬ 
pation numbers. From these results it is clear that for 
the band insulator Si it is sufficient to optimize the oc¬ 
cupation numbers to increase the band gap in line with 
experiment; the KS states are evidently already a very 
good representation of the natural orbitals. These results 
are in line with our experience with finite systems which 
shows that orbital optimization results roughly up to 25% 
of the total correlation energy and the rest 75% comes 
from the occupation numbers optimization. As may be 
seen in Fig. the opposite situation holds for the case of 
the Mott insulators NiO and MnO: clearly the KS states 
differ profoundly from the natural orbitals. In this case 
it is crucial to optimize the natural orbitals. The reason 
for this is that in the case of Mott insulators it is the 
localization of electrons which leads to formation of the 
gap and KS orbitals are not sufficiently localized, thus 
optimization over the natural orbitals is required. 

A conhrmation of this charge localization may be seen 
in the charge density. In Fig. we plot the dif¬ 
ference p(r) — plda[y), for (i) RDMFT (lower panel) 
and (ii) the LSDA-|-t7 functonal[42] (upper panel) within 
DFT for NiO. LSDA-I-C7 method is chosen because, like 
RDMFT, it also finds the correct insulating ground state 
for NiO[Sll33]. The LSDA-|-t7 method achieves this by 
both spin order an on-site Hubbard U and, in contrast to 
RDMFT, cannot treat the non-magnetic insulating state 
of this material. The impact of this difference on the 
charge density is clear in Fig. significant charge lo¬ 
calization is seen only in the RDMFT density. Inter¬ 
estingly, one observes an almost spherical charge accu¬ 
mulation at the oxygen site, a result in agreement with 
experiment [44] . but different from that found in the cor¬ 
responding LSDA-|-t7 result. 


SUMMARY 

To summarize, in this work we demonstrate that 
RDMFT in conjunction with the power functional is able 
to capture the insulating state of NiO and MnO in ab¬ 
sence of long range spin order. However, while spin order 
does not drive the insulating ground state, the large lo¬ 
cal moments in these materials require spin be explicitly 
taken into account for excellent agreement with exper¬ 
imental spectra to be obtained. The power, a, in the 
power-functional is an indicator of the amount of corre¬ 
lation and a detailed analysis shows a highly non trivial 
behaviour of the spectra, for Mott insulators, as a func¬ 
tion of a, which must lie within a small range (between 
0.65 and 0.7) for the correct insulating ground-state is 



LSDA+U 




FIG. 4: (Color online) Difference between the LSDA charge 
density and charge density calculated using LSDA-I-U and 
RDMFT, (p(r) — Plsda{jc)) for NiO. Positive values indicate 
localization of charge as compared to LSDA. 


obtained. It is further shown that the natural orbitals 
for the strongly correlated materials, NiO and MnO, are 
much more localized as compared to the Kohn-Sham or¬ 
bitals, which enables them to capture the physics of Mott 
localization in these materials. 
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